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ABSTRACT 

We investigate the relationship between Low Mass X-ray Binaries (LMXBs) and 
globular clusters (GCs) using UKIRT observations of M31 and existing Chandra, 
XMM-Newton, and ROSAT catalogues. By fitting King models to these data we have 
estimated the structural parameters and stellar collision rates of 239 of its GCs. We 
show a highly significant trend between the presence of a LMXB and the stellar colli- 
sion rate of a cluster. The stellar collision rate is found to be a stronger predictor of 
which clusters will host LMXBs than the host cluster mass. We argue that our results 
show that the stellar collision rate of the clusters is the fundamental parameter related 
to the production LMXBs. This is consistent with the formation of LMXBs through 
dynamical interactions with little direct dependence on the neutron star retention 
fraction or cluster mass. 
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1 INTRODUCTION 

Since early studies of X-ray sources in the Milky Way it 
. has been kn own that GCs are a rich source of LMXBs (e.g. 
' IClarkI I1975I ) . To explain th e populations observed in the 
Milky Way l|Liu et al.ll200ll . and references within), its GCs 
need to be two orders magnitude more efficient at forming 
LMXBs than other Galactic regions. Theorists have long 
argued that the likely reason for this is that GC LMXBs 
are formed mainly through dynamical interactions instead 
of the evolution of primordi a l binaries (e.g. Clark 1975; 
Fabian, Pringle, & Rees' 'l975|; ISutantvol Il975l : iHills fc Davl 
1976. : Verbunt fc Hut 1987 ) . The dynamical processes which 



review of lFabbianoll2006l ). From these studies it has been 
possible to identify relationships between the properties of 
a cluster and the probability that it will contain a LMXB. 

In both the Milky Way and nearby galaxies it has 
been found that LM XBs favour brighter (and hence more 



G Cs (e . g. [Sarazin et al.l l2003l: iKim et al.1 I2OO6I : 

" 3Q7D- 



can form LMXB systems include the neutron star capturing 
a donor star through tidal capture; exchange interactions; 
and direct collisions. These interactions are likely to be very 
common in the cores of GCs where the stellar densities are 
relatively high. Studies of the 13 bright LMXBs in the Milky 
Way GCs are consistent with dynamical formation scenar- 
ios. 

In the Chandra era it has been possible to study much 
greater numbers of LMXBs by looking at extra-galactic X- 
ray sources. Observations of nearby galaxies generally show 
large numbers of LMXBs associated with their GCs (e.g. the 
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ISmits et al.l I2OO6I : iKundu. Maccarone. fc Zepj I2007D . The 

likely reason for this is that higher mass clusters will 
generally have more stellar interactions and therefore 
form more LMXBs through dynamical interactions. It has 
also been suggested that higher mass clusters will re- 
tain more of the neutron stars they produce. This is be- 
cause more massive clusters will generally have higher es- 
cape velocities and neutron stars may be formed with 
large kick velocities. However it is possible for neutron 
stars to form with lower velocities via electron-capture 
supernovae (e.g. Pfahl. Rappaport. fc Podsiadlowskil I2OO2I : 



llvanova et"aLll200a ) 

If dynamical formation is the primary method of form- 
ing LMXBs in GCs, we expect there to be a strong cor- 
relation between the stellar collision rate of a cluster and 
the presence of a LMXB. The stellar collision rates in 
GCs will be dominated by collisions in their cores (where 
the stellar densities are highest). Unfortunately, due to 
their small sizes and large distances from us, resolving 
the cores of extra-galactic GCs requires very good spa- 
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tial resolution. Some previous work has tried to infer col- 
lision rates without resolv ing the cluster cores (e.g. in 
M31; iBellazzini eraP Il995h . but direct relationships be- 
tween LMXBs and collision rate have been limited to 
GCs in the Milky Way and a few nearby galaxies ob- 
served with the Huhhle Space Telescope (HST). HST ob- 
servations have been used to estimate the structural pa- 
rame t ers of GCs arou nd M31 ( Barmbv. Holland, fc Huchral 
20021: iBarmby et a,l.l [2007), Cen A ([Harris et al.l 12001 



Jordan et al 
l|jordan et al.l 



200' 



iMcLaughlin et al.l |200^ ) and M87 
20041 ) . From these studies correlations be- 



tween stellar collision rates and the presence of LMXBs have 
been proposed in Cen A and M87. 

In this letter we study 15 LMXBs in the M31 GC sys- 
tem. These are the closest extra-galactic GCs and the only 
ones for which reliable core radii can be estimated using 
ground based photometry. 



3 M31 GLOBULAR CLUSTERS 
3.1 Selection of GCs 

Throughout this letter we select GCs using the Revised 
Bologna Catalogue (RBC) of M31 GCs v3.5 and adopt the 
naming conventions of this catalogue (|Galleti et al.l |2003 . 
I2OO6, 2007). The RBC is based on the original catalogue 
of iBattistini et al.l (IiQSTI ) and is regularly revised and sup- 
plemented to include new results on the GC system. We 
identify clusters in our images based on their locations in 
this catalogue. To minimise contamination of our sample 
we only include GCs present in the RBC which are marked 
as confirmed and which were identified from the work of 
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(I992I); 


Mocheiska et al.l ( 
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(I2OOOI): 
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total 239 clusters are located in our observations with 115 
present in multiple images. 



2 DATA 

To investigate the properties of the GCs in M31, we obtained 
K band photometry using the Wide Field CAMera (WF- 
CAM) on the United Kingdom Infrared Telescope (UKIRT) 
under the service program USERV1652. These data cover 
most of the M31 disk, but avoid the central regions where 
surface brightness fluctuations are largest. The data include 
239 of the confirmed M31 GCs. 

Each observation consisted of microstepped 5s expo- 
sures with a nine point jitter pattern. The 2x2 microstep- 
ping was used to give an effective pixel size of 0.2" and ensure 
the images were well sampled. Five observations were taken 
of each field to give a total exposure time of 225s. This gives 
a detection limit of K=19 at 5(t. 

The images were processed using the standard WF- 
CAM p ipeline (for detai ls on the WFCAM and its pipeline 
see e.g. iDve et aT]|2006l ). The pipeline reduces and stacks 
the raw images and adds an accurate astrometric solution 
to the final images by fitting sources to the Two Micron 
All Sky Survey (2MASS). For combining our microstepped 
images we chose not to use the standard pipelines inter- 
leaving method. Instead th e reduced images were drizzled 
together (|Fruchter fc Hookl l2002) using the IRAF STSDAS 
task DRIZZLE. Drizzling has the advantage that it produces 
combined images with smoother PSFs than interlacing. 

To measure the structural parameters of the GCs ob- 
served, spatial resolution is very important. Our observa- 
tions were taken on the nights of 2005-11-30 and 2007-08-06 
with seeing of 0.85-1.00" and 0.6-0.8" respectively. At the 
distance of M31, this corresponds to a spatial resolution of 
3.2-3.9pc and 2.1-2.9pc. This is better than can be achieved 
by HST observations of Virgo Cluster galaxies where the 
spatial resolution is ~8pc and only slightly worse than Cen A 
where the resolution is ~1.8pc. Our observations also have 
the advantage over these HST observations that we are able 
to measure our PSF directly using stars in our fields instead 
of having to rely on models. 



3.2 Profile fitting 

To investigate the structure of the GCs in our observations, 
we must first consider the effects of seeing on their appear- 
ance. The images of the GCs studied are the result of a 
convolution of their physical size with the PSF of the obser- 
vations. To account for this, we create a model for the PSF 
of each image so that we can deconvolve it from the clusters 
and investigate their underlying properties. The PSF was 
modelled with a Moffat profile using the DAOPHOT tasks 
ALLSTAR/SEEPSF to fit the brightest 40 stars which were: 
unsaturated; had no bright neighbours; and were greater 
than 100 pixels from the detector edges (where the noise is 
significantly higher) . No significant variation in the PSF was 
observed across the images so we select a single PSF model 
for each detector of each observation. 

To find the structural param eters of th e GCs we fit 
spherical, single mass King models (iKindl 19661 ) to their pro- 
files using the program SUPERKING by Waters et al. (in 
prep.). The full details of this fitting and analysis will be pre- 
sented in a subsequent paper. Here we present the results 
of the fits along with some consistency checks to test the 
reliability of the parameters found. The models were first 
convolved with the appropriate PSF for the image before 
being fit to the cluster profile. The best fitting model was 
then selected based on minimisation and the structural 
parameters for this model were output. 

Before fitting the GC profiles, bright stars (above 5(7) 
were removed from the GC region using the DAOPHOT task 
ALLSTAR. Each cluster was fit out to a radius of 20" . This 
was chosen to extend beyond the expected tidal radius for 
85% of the clusters (based on Milky Way GCs in the Harris 
catalogue; Harris 1996). This allows for the tidal radius and 
background level to be computed as accurately as possible. 



3.3 Derived parameters 

Table [1] shows the best fitting structural parameters for the 
239 GCs studied. Where clusters were present and well fit in 
more than one observation, we select the observation with 
the best seeing. We assign a fiag to each cluster based on 
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Table 1. Catalogue of M31 GC parameters 
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Table is available in full as Supplementary Material in the ele ctronic edition of the journal or from the VizieR archive 

1 Names taken from the RBC of M31 GCs l lGalleti et al.ll2007^ 

2 Assuming the distance of M31 to be 7 80kpc i McConnachie et al.ll2005h 



X-ray source associated with GC from; lTrudolvuboT^r?riedhorsk^^004h fT04'l: ISuDPe7 et al ] l|2001^ (SOI); or no known source (0) 



3j 

^Reliability flag based on visual examination of the cluster and its best fltting model. 



visual examination of the cluster and its residual after sub- 
traction of our best fitting model. These flags are used to 
identify poorly fitting models, as well as clusters contami- 
nated by very bright stars. 

Using these parameters we estimate the stellar collision 



rate (F) for each cluster via T (x p, 



3/2^2 



This relationship 



assumes a constant collision rate per unit volume in the 
cluster cores which dominates the total collision rate. It fol- 
lows from the more general formula F oc p^r^/a by assum- 
ing Virial e quilibrium to estimate the velocity dispersion as 



a oc p'J'^rc jVerbunt &: Hutlll987t ). This relationship is used 
to calculate the collision rates presented in table [T] 

Despite studying only the confirmed clusters in the 
RBC it is possible that our sample contains some contamina- 
tion from stellar sources. We identify misclassified clusters 
by finding objects whose stars should have collided many 
times over a Hubble time. We can estimate the timescale for 
stellar collisions from F by assuming the cluster is comprised 
of solar type stars in orde r to find the constant in the above 
proportionality (see eq. S iVerbunt fc Hutlll987^ . This gives 
that the timescale on which we expect a star to collide is: 



rv 



1.1 X 10"* 



pc- 



105yr-i 



(1) 



Ft is the total collision rate and Ni, = {A/2>)-Kr^pc is the 
total number of stars in the cluster core. From this equation 
we find 5 GCs to have tcoii < lO^yr. One of the clusters 
(B431) has a bright star contaminating its profile. However 
the other four clusters (B006D, B096D, B292D, and B480) 
are well fit and we propose that they are likely to be mis- 
classified stars. It has previously been suggested that B292D 
may not be a GC (|Huxor et al.| [2008) and we note that the 
radial velocity estimates in the RBC of three of these clus- 
ters are consistent with being Milky Way halo stars. We 
remove these clusters from the following analysis. 



3.4 Consistency checks 

Having identified the best fitting King model to describe 
each cluster we investigate the reliability of the parameters 
found. Due to the small spatial scale of the clusters we are 
fitting, it is likely that the errors on these parameters will be 
dominated by errors in the PSF model which is convolved 
with the cluster. 

To investigate the magnitude of the typical measure- 
ment errors, we compare the results obtained from fitting 
the same cluster from different observations. In total, 115 



10 



-0.5 



+ ^ 



++ + + 



" , ^■■■■^iT^--+-**ff*ri*«t^--^ 



16 15 14 13 12 

K (mag) 

Figure 1. Comparison of the parameters derived from fltting the 
same cluster in different observations. 



clusters are present in multiple images. These observations 
have all been taken under slightly different conditions and 
will have different PSF models. Therefore by fitting these 
clusters independently and comparing the resulting param- 
eters we can estimate the reliability of the parameters cal- 
culated. Figure [T] shows the differences between the derived 
parameters for the clusters fit in more than one image. 

For clusters brighter than K=15mag we find good agree- 
ment between the parameters found. However it can be 
seen that the scatter increases significantly for clusters with 
K>15mag. This suggests that the errors on the parameters 
found for these faint clusters are significantly higher. The 
large deviation in the magnitudes found for these faint clus- 
ters is likely to be due to contamination from nearby stars 
or surface brightness fiuctuations. 

We expect that the variation observed in figure [T] gives 
a reasonable estimate of the errors on the parameters found 
for all the GCs studied. In the following analysis, we include 
all 235 GCs studied but note that the errors on individual 
parameters for the faintest clusters will be large. Since most 
of the GCs studied and all X-ray GCs are brighter than 
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Figure 2. Comparis on of the parameters measured liere witti 
those published by iBarmbv et"al] l l2007t) (B07) using i?ST 
STIS/ACS (solid circles) and WFPC2 (open circles) observations. 



this, our conclusions should not be sensitive to this increased 
error. 
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Figure 3. Stellar collision rate vs magnitude for all GCs (crosses) 
and LMXB hosting clusters (solid diamonds). This shows a clear 
relationship between the two parameters which is in good agree- 
ment with the predicted relationship F oc M^^^ (dashed line). 
The LMXB hosting clusters are found to have higher than aver- 
age collision rates for their magnitude. 



a clear scaling difference. This difference is to be expected 
and is primarily because we estimate the number of stars 
based on K-band instead of V-band photometry. This off- 
set is therefore most likely due to the colour of the clusters. 
Since our analysis is relative only to values measured in this 
study, this offset should not effect our conclusions. Also since 
stellar mass is more closely correlated with K-band luminos- 
ity than any optical band, comparisons of the collision rates 
within our data are likely to be robust. 



3.5 Comparison with previous work 

Some of the GCs in M31 have been observed by the 
HST under several different programs. Using these data 
the structural parameters for 96 clusters have already 
been esti mated from either ACS, STIS or WFPC2 obser- 
vatio ns (|Barmbv. Holland, fc Huchral |2002| : iBarmbv et al.l 
Figure [5] compares these previously found parame- 
ters with those found in this study. There are 33 clusters 
present in both datasets. 

For clusters brighter than K=15mag, we find reasonable 
agreeme nt between the radii found in this study and those 
found bv lBarmbv et all (|2007l ). For clusters fainter than this 
we find significant differences between both the core and 
tidal radii found. However this is within the large errors 
predicted in the previous section for these fainter clusters. 

For all the cluster parameters compared we find greater 
discrepancies between our results and those based on 
WFPC2 observations (open circles). The likely reason for 
this is the methods used to estimate these parameters. For 
our data, we fit the whole profile of the cluster and select 
the best fitting tidal radius. While the WFPC2 observations 
are fit to the inner region of the cluster (where the signal 
to noise is greater ) and the tidal radius i nferred from the 
best fitting model (IBarmbv. HoUand. fc Hu chra 2002). This 
second method is much more susceptible to slight deviations 
from a pure King model. We therefore believe our tidal radii 
may better represent the actual tidal radii of the clusters. 

The stellar collision rates found by both studies show 



4 X-RAY SOURCES IN M31 GCS 

X-ray sources have bee n associated with G Cs in M31 from 
ROSAT observations (jSupper et al.l 200ll) and archived 
Chan dra and XMM-Newton i Trudolvubov fc Priedhorskvl 
|2004 ) observations. For details on the depth, coverage and 
reliability of these catalogues we refer the reader to these 
papers. Of the GCs currently known to host LMXBs, 15 are 
included in our data. 



4.1 LMXB relationships 

Figure[3]shows the stellar collision rate as a function of mag- 
nitude for all GCs studied (crosses) and those containing a 
LMXB (diamonds). It can be seen that the LMXBs favour 
brighter clusters and those with higher stellar collision rates. 
A Kolmogorov-Smirnov (K-S) test between all GCs and the 
LMXB hosting GCs shows that both of these relationships 
are significant with probabilities of 10^^ for luminosity and 
10^^ for collision rate that they come from the same popula- 
tions. This demonstrates that the collision rates and masses 
of GCs are very good discriminators in selecting LMXBs. 

It can be seen from figure|3]that there is a clear relation- 
ship between the mass of a GC and its stellar collision rate. 
In order to explore the relative effects of these parameters 
on LMXBs, we must first consider this relationship. Assum- 
ing the relationship to be linear between the logarithm of 
these values, we find that F cx Ml^f. This is consistent with 
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Cumulative Fractional Collision Rate 

Figure 4. This shows the cumulative fractional collision rate for 
all GCs (small crosses) and LMXB hosting GCs (large crosses). If 
the formation of LMXBs is linearly proportional to GC collision 
rate, they should be evenly distributed along this plot. 

the Milky Way GCs and the theo retical approximation that 
r oc Mtl (|Davies. Piotto. fc de A ngcli 2001). 

It can be seen that the LMXB hosting GCs have higher 
than average collision rates for their mass. To investigate 
whether this is a statistically significant effect, we first de- 
trend the data using the relationship between F and mass. 
We then run a K-S test between all GCs and the LMXB GCs 
and find a probability of 10~^ that they are from the same 
distribution for a detrending on F oc Ml^^. We note that 
there is some uncertainty in the actual relationship between 
F and mass due to the large scatter in figure [S] To ensure 
our results are robust to errors in this relationship, we re- 
run the tests assuming F oc M/^f^'^'^^ and find probabilities 
of 10~^ and 10~^ respectively that they are from the same 
distribution. This demonstrates that, even for the steepest 
reasonable relationship, the LMXBs favour GCs with higher 
than average stellar collision rates for their mass. 

To investigate whether the formation of LMXBs is lin- 
early proportional to t he stellar collision rate of a cluster, 
we follow the method of IVerbunt fc HutI l| 19871 ). First we di- 
vide all GC collision rates by the total collision rate for the 
whole GC system. We then sort all the clusters by this frac- 
tional collision rate and find the cumulative value for each 
cluster, such that the cluster values run from to 1. The 
result of this is that the total collisions occurring in the GC 
system should now be evenly distributed between and 1 
(i.e. we expect 10% of the total collisions to occur in clusters 
with values in the range 0-0.1). Therefore, if the formation of 
LMXBs is linearly proportional to collision rate, their host 
clusters should be evenly spaced in this plot. It can be seen 
from figure|4]that the LMXB hosting clusters (large crosses) 
are consistent with this. 



5 CONCLUSIONS 

We show clear relationships between the presence of a 
LMXB and both the stellar collision rate and mass of its 
host GC. 

The stellar collision rate is found to be the best discrim- 
inator in selecting LMXB hosting GCs. We suggest that the 
weaker relationship between mass and LMXB presence may 
be primarily due to the relationship between mass and stel- 
lar collision rate. Our results demonstrate that the stellar 
collision rate is likely to be a fundamental parameter related 
to the formation of LMXBs. This result is in agreement with 
previous studies of GCs in the Milky Way and Cen A. 

The linear relationship found between the presence of 
LMXBs and stellar collision rate is in good agreement with 



the systems being formed by dynamical interactions. This 
also suggests that the current dynamical properties of the 
GCs are related to their current LMXB populations. 
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